Exciting pions in nuclei: DCX and Electroproduction in the resonance region 
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We describe the double charge exchange (DCX) reaction of pions with different nuclear targets 
at incident kinetic energies of 120 — 180 McV within the GiBUU transport framework. The 
DCX process is highly sensitive to details of the interactions of pions with the nuclear medium 
and, therefore, represents a major benchmark for any model of pion scattering off nuclei at low 
and intermediate energies. We achieve a good quantitative agreement with the extensive data set 
obtained at LAMPF. Furthermore, we present first results on electron scattering off nuclei for beam 
energies of 700-1500 MeV and virtualities of Q 2 < 700 McV 2 . Including quasi-elastic and pion- 
production processes and employing in-medium kinematics and realistic form factors, we evaluate 
inclusive cross sections and compare to data obtained at the ADONE storage ring (Frascati). 
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I. INTRODUCTION AND MOTIVATION 

Major experimental efforts arc being undertaken to un- 
derstand the in-medium properties of hadrons, e.g. |l|; 
0; [||). In many cases, the hadrons of interest decay in- 
side the nuclear medium while one finally observes its 
reaction products outside the medium. However, one 
must discriminate between the change of the original 
hadron properties in the medium and many-body final 
state re-scattering effects which influence the observed 
decay products in the case of the hadronic observables. 
In this context, transport models become relevant: these 
models describe the whole reaction process in a micro- 
scopical manner by evolving the phase space density of 
all involved particle species in time. Therefore, the trans- 
port approach is suited to shed light on this discrimina- 
tion problem. 

The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) 
transport model j3) has been developed over the last 20 
years to describe heavy ion collisions, photon-, electron-, 
pion- and neutrino-induced reactions within one unified 
transport framework. After a major restructuring phase 
of the whole source code, we report here new applications 
of this method. 

In this work we focus on low-energy processes, where 
pions play a dominant role. Recently, we have shown 
using GiBUU (@; H) that pion re-scattering in the final 
state description of photon induced double-pion produc- 
tion produces a considerable modification of the itit in- 
variant mass distributions observed by the TAPS collab- 



oration (P). The propagation of low-energy pions in nu- 
clear matter in the GiBUU framework has already been 
extensively discussed in ref. (@) and compared to quan- 
tum mechanical calculations. In the spirit of benchmark- 
ing the model, pionic double charge exchange (DCX) is 
another very interesting reaction. The fact that DCX re- 
quires at least two nucleons to take place makes it a very 
sensitive test for pion re-scattering and absorption. 

We also present our first results for electron scatter- 
ing off nuclei for beam energies between 700 MeV and 
1500 MeV and virtualities of Q 2 < 632 MeV 2 . Employ- 
ing impulse approximation, we obtain inclusive cross sec- 
tions with full in-medium kinematics and realistic form 
factors for quasi-elastic scattering and pion production. 

This article is structured in the following way. First 
we introduce our GiBUU transport model emphasizing 
the most relevant issues. Next, we discuss two distinct 
applications: DCX and electron scattering off nuclei. 



II. THE GIBUU TRANSPORT MODEL 

Boltzmann-Uehling-Uhlenbeck (BUU) transport mod- 
els are based on the Boltzmann equation, which was mod- 
ified by Nordheim, Uehling and Uhlenbcck to incorporate 
quantum statistics. A brief description of the formal- 
ism is given below. For a detailed discussion concerning 



the physical input for pion-induccd reactions we refer the 
reader to f(| 0) < 



and the references therein. 



* Elect ronic address: [oliver.bussOthco, physik. uni-gie ssen.de| 

URL: |http://th eorie . physik . uni-giessen . de/~ ol iver| Work 
supported by DFG. 

tAlso in RRC "I.V. Kurchatov Institute", 123182 Moscow, Russia. 
Present address: Frankfurt Institute for Advanced Studies (FIAS), 
Germany. 



A. The BUU equation 

The BUU equation actually consists of a series of cou- 
pled differential equations, which describe the time evolu- 
tion of the single-particle phase-space densities f a (r,p, t). 
The index a = tt, w, N, A, . . . denotes the different parti- 
cle species in our model. A large number of mesonic and 
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baryonic states is actually included, but at the energies 
of interest for this study, the relevant ones are tt, N and 
the A(1232) resonance. 

For a particle of species X, its time evolution is given 
by 

dfx (r, P, t) dfx (r, p, t) dH x dfx (f, p, t) 



dt dt 
dH x df x (r,p,t) 



dr dp 
with the one-body Hamilton function 



dp dr 

= Icoll(fx,f a ,fb,...) (2.1) 



Comparing to a momentum independent nuclcon poten- 
tial, which is approximately —50 MeV deep, the A po- 
tential is, therefore, taken to be 



A° a (p,7) 



^nucleoli (P> 



Here we assume the same momentum dependence for the 
nucleon and the A potential. No potential for the pion 
was included in the calculations. 



III. PIONIC DOUBLE-CHARGE-EXCHANGE IN NUCLEI 



Hx(r,p) 



p - A x (r,p)) + m\ + U x (r, p) 



+A° x ^p) 



(2.2) 



The scalar potential Ux and the vector potential A x of 
species X may in principle depend upon the phase space 
densities of all other species. Hence, the differential equa- 
tions are already coupled through the mean fields. In 
the limit of I co u = 0, cq. (|2.ip becomes the well-known 
Vlasov equation. The collision term I co u on the right- 
hand side of eq. (|2.1[) incorporates explicitly all scatter- 
ing processes among the particles. The reaction proba- 
bilities used in this collision term are chosen to match 
the elementary collisions among the particles in vacuum. 
Therefore, we model these probabilities as a sum of dom- 
inant resonance terms and small background terms. The 
latter resonance terms are based on the partial wave anal- 
ysis of Manley et al. (d) . The background terms are de- 
fined as differences between experimental data and reso- 
nance contributions. Note that background contributions 
are instantaneous in space-time, whereas the resonances 
propagate along their classical trajectories until they de- 
cay or interact with one or two nucleons in the medium. 

Concerning pion absorption, the most important mech- 
anisms are a two step process in which irN — > A is fol- 
lowed by AiV -> NN or ANN -> NNN, and the one 
step background process irNN — > NN. 



B. Medium effects 

As medium modifications, we include Pauli-blocking 
and Fermi- motion of the nucleons and Coulomb forces. 
Furthermore, we account for the collisional broadening 
in the A resonance width. The hadronic potentials in- 
cluded are introduced as time-like components of vector 
potentials in the local rest-frame (0). For our purposes 
the most important mean field potentials are those act- 
ing on the nucleon and the A resonance. The nucleon 
mean-field potential is parametrized by Welke et al. (fiot ) 
as a sum of a Skyrme term depending only on density and 
a momentum-dependent contribution (for explicit details 
and parameters see (0)). Phenomenology tells us that the 
A potential has a depth of about —30 MeV at po (fill ; [TT 



A. Introduction 

The double charge exchange (DCX) reaction of pions 
received a considerable attention in the past (see for in- 
stance Ref. (p~3h and references therein). The mechanism 
of two sequential single charge exchanges has tradition- 
ally been able to explain the main features of this reac- 
tion at low energies although the contribution of 
the A(n, ttit)X reaction becomes progressively important 
as the energy increases (fill [P7T ). At higher (~ 1 GeV) 
energies, the sequential mechanism becomes insufficient 
to account for the reaction cross section (fl8l ; fl9h . Exten- 
sive experimental studies performed at LAMPF obtained 
high precision data for doubly differential cross sections 
on 3 He (flO) and heavier nuclei ( le O, 40 Ca, 208 Pb) Q 
in the region of E kin = 120 - 270 MeV. 

Huffier and Thies ((22] ) explored for the first time the 
applicability of the Boltzmann equation in irA collisions 
and achieved qualitative agreement with data on single 
and double charge exchange. Their method to solve the 
Boltzmann equation was based upon an expansion of the 
pion one-body distribution function in the number of col- 
lisions. There, in contrast to our work, the Boltzmann 
equation is not solved with a test-particle ansatz but by 
reformulating it into a set of coupled differential equa- 
tions which can then be solved in an iterative manner. 
However, this approach was based on simplifying assump- 
tions of averaged cross sections and averaged potentials. 
The work by Vicente et al. (fl6r l was based upon the cas- 
cade model described in (f23h. There, a microscopic model 
for ttN scattering was used as input for t he p ion reaction 
rates in the simulation. In that work (fig ), pion DCX 
off 16 O and 40 Ca has been explored and fair quantitative 
agreement with data has been achieved. 

In the following, we explore DCX on heavier nuclei, 
comparing with the data measured by Wood et al. |2ll). 
We also address the scaling of the total cross section dis- 
cussed by Gram et al. ((H). To focus only on single-pion 
re-scattering, we consider incoming pion energies below 
^kin = 180 MeV; above that energy 2ir production be- 
comes prominent and DCX does not happen necessarily 
in a two-step process anymore. Due to the small mean 
free path of the incoming pions, the process is mostly 
sensitive to the surface of the nucleus. In (0) we also ad- 
dressed the issue of parallel versus full ensemble scheme 
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in the numerical implementation of the BUU equation. 
There, we found that the parallel ensemble scheme is a 
good approximation and has therefore been applied in 
the following calculations. 



B. Influence of the density profile 

Due to the low pion mean free path in nuclear mat- 
ter (@), DCX is very sensitive to the surface properties 
of the nuclei. Therefore, we compared the results with a 
neutron skin for the Pb nucleus to the results obtained 
without such a neutron skin; for details see (0). Neutron 
skins are very interesting because in those skins only ir + 
mesons can undergo charge exchange reactions. For the 
positive pions this causes an enhancement of DCX pro- 
cesses at the surface, so the pions do not need to pen- 
etrate deeply for this reaction. Hence, the probability 
for their absorption is reduced. As can be observed in 
fig. 1, the enhancement in the total cross section for 
Rn - R P = 0.266 fro. is roughly 35% at 180 MeV. The 
accurate determination of neutron skins is relevant for 
different areas of physics such as nuclear structure, neu- 
tron star properties, atomic parity violation (PV) and 



% A -> X ji + 



heavy ion collisions (|25l : 126 
periment (PREX) at JLab 



J27|). The Parity Radius Ex- 
28j) shall measure the neutron 
radius with high precision (1%) using PV electron scat- 
tering. We can see that the DCX cross section is very 
sensitive to the size of the difference between the proton 
and the neutron radii of Pb. The effect is specially large 
(more than a factor 2) at forward angles, where our model 
performs very well (see next Section). Indeed, without 
neutron skin, due to strong pion absorption in the bulk of 
nuclear matter, the DCX cross section is small at forward 
angles. The presence of a neutron skin favors DCX in 



n + Pb -> X n 
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FIG. 1 Influence of the density distribution on the angular 
distributions for the double charge exchange process tt + Pb — > 
7T~ X at _Ekin = 180 MeV. The solid line shows the result 
where we implement different radii of the density distribu- 
tions of neutrons and protons in our model (Rn — R v — 0.266 
fm). As a consequence, the nucleus is surrounded by a small 
neutron skin. For the dashed line we assumed same radii for 
the proton and neutron distributions. 
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FIG. 2 The inclusive double charge exchange total cross sec- 
tion as function of the nuclear target mass at E^ n = 120, 150 
and 180 MeV. The lines connecting our results are meant to 
guide the eye; the data are taken from Ref. (plT ) (lower panel: 
Eun = 120, 150 and 180 MeV, upper panel: only 180 MeV). 



peripheral reactions, where the pion propagates in prac- 
tically pure neutron matter. This naturally enhances the 
DCX cross section at forward angles. Hence, a precise 
measurement of DCX at forward angles, combined with 
a realistic theoretical analysis could be a valuable source 
of information on the neutron skins complementary to 
the one obtained with PV electron scattering. Note that 
for the 7T~ similar arguments lead to a reduction of the 
cross section. 



C. Comparison to data 

Now we proceed to the comparison with the data mea- 
sured at LAMPF by Wood et al. Q. We discuss first 
the total cross section and thereafter the angular distri- 
butions. 1 

In fig. [2] one can see the excellent quantitative agree- 
ment to the total cross section data at 120, 150 and 
180 MeV for Oxygen and Calcium. For the Lead nucleus 
we see some discrepancies, which are, however, within 



In (0) we also show double differential cross sections as a function 
of both angles and energies of the outgoing pions. 
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FIG. 3 The scaling of the total charge exchange cross section 
according to eq. (|3.1[) is visualized by dividing atot by the fac- 
tor A 2 / 3 Q(Q~1)/(A-1) and plotting it as afunction of A-Q. 
Q denotes the number of protons in the case of (7T - , tv + ) and 
the number of neutrons in (7r + ,7r _ ). The points are GiBUU 
results at pion kinetic energies of 180 MeV; the dashed line 
denotes a function proportional 1/(A — Q), corresponding to 
the exact scaling. 



the experimental accuracy. Notice that we reproduce the 
different A dependencies of both (tt + 7 it~) and (n~,n + ) 
reactions. It is due to the fact that, when A increases, the 
number of neutrons increases with respect to the number 
of protons, and this favors the tt + induced reaction. 

In Hil), Gram et al. discuss a scaling law of the to- 
tal cross section. They argue as follows. Since the first 
collision takes place predominantly at the surface, the 
cross section should scale with A 2 / 3 . Furthermore they 
assume that DCX is mainly a two-step process, and that 
a pion which undergoes an elastic process at the first 
collision will not contribute. This is reasonable because 
the incoming pions loose energy in the elastic process, 
and their cross section for a second charge-exchange re- 
action is hereafter very much reduced. For a negative 
pion the first charge exchange reaction occurs with a 
probability of Z/N where Z(N) denotes the number of 
protons (neutrons). This is the case if the interaction is 
dominated by the A resonance, as it should be in this 
energy region. Finally, the second charge exchange pro- 
cess then takes place with the probability (Z—1)/(A—1) 
since, in the isospin limit, the ir° interacts equally well 
with protons and neutrons. Putting these considerations 
together and extending them to the ir + case, the cross 
section for DCX is expected to scale according to 



O"tot 



A 2/3_ 



Q Q-l 
A-QA-1 



(3.1) 



where Q denotes the number of protons in the case of 
7r~ induced and the number of neutrons in n + induced 
DCX. 

Gram et al. (24) find good agreement of this scal- 
ing law with experimental data. This scaling is fulfilled 
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FIG. 4 Angular distributions for the double charge exchange 
process tt+A -v ty^X at _B kin = 120, 150 and 180 MeV. 
The data points are taken from |2ll) ; only systematical er- 
rors are shown. The solid lines represents the GiBUU results 
at 120 MeV, the dashed lines at 150 MeV and the dotted lines 
at 180 MeV. The bold lines represent the results while the 
thin lines represent a la confidence level for each point based 
upon our statistics. 



within our GiBUU simulations as can be seen in fig. [3] 
Nevertheless, one may wonder why this scaling law works 
in a process which is so sensitive to the neutron skin on 
heavy nuclei, as has been shown in fig. 1. Since the 
first collision takes place on the surface, a neutron skin 
causes an enhancement in the A(tt + , it~)X reaction while 
A(ir~ ,ir + )X is suppressed. This effect leads to a devia- 
tion from the scaling. However there are also Coulomb 
forces which arc not negligible. The Coulomb force en- 
hances A(ir~ , tt + )X by attracting the negative projectiles 
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FIG. 6 Same as figHfor 208 Pb. 



good quantitative agreement for both O and Ca. In the 
Pb case, the (7r _ ,7r + ) reaction is well described, but the 
(7r + ,7r~) one is underestimated in spite of the enhance- 
ment caused by the neutron skin. 



FIG. 5 Same as figHfor i0 Ca. 



and repelling the positive products, which therefore have 
a smaller path in the nucleus and undergo less absorption. 
And, due to similar arguments, the reaction A(ir + , tt~)X 
is suppressed. We find that this effect counteracts the 
one from the neutron skin restoring the scaling. In any 
case, the approximate scaling exhibited by the cross sec- 
tion shows that the reaction is very much surface driven 
and can be very well understood in terms of a two-step 
process. 

In fig. HGUwe show da/dVL for DCX at E kin = 120, 150 
and 180 McV on le O, 40 Ca and 208 Pb as a function of the 
scattering angle 9 in the laboratory frame. Our results 
(bold lines) are shown together with their uncertainties 
of statistical nature (thin lines) . The latter ones are well 
under control except at very small and very large angles, 
where statistics is very scarce. Again, there is a very 



IV. ELECTRON INDUCED PROCESSES 

Our investigation of electron induced processes in the 
resonance regime has two major goals. On the one hand 
side, we want to investigate the influence of final state 
events on possible electron signals; and on the other hand 
side we want to understand the role of the hadronic res- 
onances in the scattering of electrons and nuclei. 

Furthermore, Leitner et al. ((291 ; l30h have already per- 
formed calculations for neutrino-induced reactions in the 
resonance energy regime using the GiBUU transport 
model. The electron and neutrino interaction with nu- 
clei is closely interconnected, so we can test assumptions 
going into the neutrino calculations. This is interesting, 
since unlike in the electron sector there are almost no 
experimental data available in the neutrino sector. 

The wavelength of the photon is considered small 
as compared to intra-nucleon distances in the nucleus. 
Hence we consider the nucleus as a Fermi-Gas of nucle- 
ons, and the total reaction rate is given by an incoherent 
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FIG. 7 (Quasi-)Elastic electron-nucleon scattering. 
l' 
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FIG. 8 Electron induced pion production off a nucleon. 



sum over all nucleons. These nucleons are under the influ- 
ence of a mean field potential as described in lll.Bl there- 
fore the dispersion relation is modified in the medium. 
However, we neglect at the present stage modifications 
of the nucleon width due to short range correlations. 



A. (Quasi-) Elastic Scattering off a single nucleon 

In this section we will treat the most trivial case of 
electron-nucleon scattering: the elastic scattering. In the 
case of electron scattering off a bound nucleon this pro- 
cess is referred as quasi-elastic (QE). 

Our notation is shown in fig. [7] the in- and outgoing 
lepton momenta are denoted I and I', the nucleon mo- 
menta p and p', the photon momentum by q = I — V . 
The 7 matrices and the spinor normalization are taken 
according to Bjorken and Drell (f3lh . The electron mass 
is denoted by m e , the masses of the in- and outgoing 
nucleons by to and to'. 

Parametrising the hadronic vertex by the operator Jq E , 
the matrix clement for this process, summed over final 
and averaged over initial spin, is given by 



T H v ^ 



with Si, Sf denoting the initial and final nucleon spin; ti, 
tf denoting the electron spins. In the latter equation, the 
electron and nucleon contributions are separated in the 
so-called hadronic and leptonic tensors Hq^ and £,,„. 



The most general form of the hadronic vertex can be 
expressed as 



QE 



2toq 



F 2 u^q v + F 3 q^ 



with the nucleon vacuum mass too ~ 0.938 GeV. De- 
manding charge conservation at the vertex, i.e. q^j^ = 0, 
we get 



^QE — e 



2to,q 



1v 



Fi(m — m') 



</" (4.1) 



where the last summand is absent in the case of non- 
momcntum-dcpcndcnt nucleon masses. In the vacuum, 
the masses of the in- and outgoing nucleon are identical, 
therefore this latter term vanishes. In the medium, how- 
ever, the masses depend through the mean fields on the 
momenta of the nucleons. The form factors i*\ and F 2 
are the standard Dirac and Pauli form factors. As an 
approximation, we choose the in-medium form-factors to 
be the same as the vacuum form factors. 

Finally, the form factors F\ and F 2 can be expressed 
by the so called Sachs form-factors G e and G m , which 
have in the electron-nucleon CM-frame an interpretation 
as charge form factor and magnetic form factor. For 
the Sachs form- factors G\ , G v ml GfL we will use the so 
called BBA03 parametrization (|32l ). For G™ we use a 
parametrization obtained by Krutov and Troitsky |33l). 
Finally, the cross section for quasi elastic scattering is 
given by the general formula 



do~ c 



1 



-drii> 



d\l'\ S((l + p - l') 2 - to' 2 ) \M\ 2 



(4.2) 



Note that we use the full in-medium kinematics, therefore 
in 7^ to' in the general case. 



B. Electron induced single-pion production off a single 
nucleon 

As depicted in fig. [51 we now consider an additional 
pion with momentum k in the final state. As in the QE- 
scattering, the matrix element can be written in the form 



\M\ 2 = L^H? 



(4.3) 



where is now including another hadronic vertex con- 
tribution as compared to the case of QE-scattcring: 

x (u{p',s f ){J^u{p,s i )) ji . 



In the same spirit as Bercnds et al. (1341 ) but in the nota- 
tion of MAID (f35h , we parametrize this hadronic vertex 
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by 

6 

J% = Y,AiM? (4.4) 

i=l 

with 

Mi = y7Vtf-*7") (4-5) 

M£ = 2z 7 5 (P^-(fc-|)-P. 9 (fc-|r) 

M£ = -i^^k-q-fW) 

M% = -2i 7 5 (7 M (Z • P - $P>*) - 2mM? 

Mq = -i^W-tf-f) 

and P = (p' + p)/2 . The form factors A\, . . . ,Aq, the 
so called invariant amplitudes, are functions depending 
on all possible scalars which one can construct in the 
vacuum out of the available 4- vectors at the vertex. As a 



minimal set of scalars we choose W — y/s, Q 2 = —q^q^ 
and the CM scattering angle 8 between q and k. The 
Ai are parametrized by the MAID group ff36k H3; l38T j 
according to elementary scattering processes. As a first 
approximation we will use the same paramctrizations also 
in the medium case. However, we have to be careful since 
W can get modified in the medium. So we define 

P v — = {sjml+f,p) 
with mo = 0.938 GeV and 



W vacuum = y((/ + p vacuum )A'(g+p vacuum ) AI . 

We choose the CM-frame of the hadronic vertex as our 
reference frame for its evaluation. Finally, we define the 
in-medium form-factors as the vacuum form-factors of 
the MAID analysis evaluated at the kinematical point 
Q 2 , W" vacuum ) The cross section in the medium is given 
by 



da* 



di ,0 dn r dn,. 



1 m 2 e m' m \k\ 2 l'° L^Hg" 

\v e -V n \ 2 (2tt) 5 fc° P '° 1° p" D 



\k\—Xn 



with 



D 



JM 
k° 



\k\ (q + p 3 )k/\k\) [l + Ari V(\p'\) + 2^p> + m 2 ^p 



V(\p'\) 



dV(\p'\) 
d\p' j 



and xo is the solution for \k\ to the energy conservation particles are transported out of the nucleus. This may, 
condition owing to re-scattcring effects, lead to a change in the final 

state particle multiplicities and distributions. However, 
q° + p° = k + p'° we will restrict ourselves in this work to inclusive cross 

~~2 sections and discuss the FSI effects in a forthcoming pub- 
lication. 

Fixing the electron kinematics, the inclusive cross sec- 
tion for scattering off a nucleus is given by 



\k\ 2 + m 2 . + \ \p + q — k\ 2 + m'(p + q — k) 



Note that m and m' are the in-medium masses and V 
denotes the potential of the nucleon, i.e. 



d(Tt, 



p' +m 2 + V(p'),p' 



C. Scattering off nuclei 

We assume that the electron scattering off a nucleus 
takes place in a two-step process. First the electron ex- 
cites a pion and a nucleon inside the nucleus, this is the so 
called initial state. In the so called final state process the 



drdQu 



J =1 



drdfiu 



PpB(rj,p') 



o'=l+p-V 




di ,0 dn v dn k 



PpB(rj,p') 



p 1 =l-\-p—V — k 



where we sum over all nucleons in the nucleus. Note also, 
that we have to integrate in the latter equation over the 
angular distribution of the pion in the — > irN con- 
tribution. The functions Ppb denotes the Pauli blocking 
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factor at initialization time. Up to now. we do not con- 
sider any other contributions than quasi-elastic scattering 
and single ir production, i.e. irir production is omitted. 
Henceforth, we expect our model to fail above the tttt 
threshold. 



D. Results 
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FIG. 9 Electron scattering off ie O at incoming electron en- 
ergies of 700 MeV and 880 MeV and an electron scattering 
angle of 37°. The plot shows a comparison of our results to 
the experimental data obtained by Anghinolfi et al. (|39[) in 
lab coordinates as a function of the virtual photon energy 
q° . The dashed curves denote the result without a mean field 
potential for the nucleon and the solid ones with such a po- 
tential. The 7*Af — > ttN contribution to each result is shown 
with additional symbols. 

In fig. [5] and [TO] we show a comparison of our inclu- 
sive cross section to the experimental data obtained by 
Anghinolfi et al. (JH). At different beam energies ranging 
from 700 to 1500 MeV, we show calculations both with- 
out (dashed curves) and with (solid curves) a mean field 
potential for the nucleon. The 7* A — > irN contribution 
to each result is shown with additional symbols. 

Especially at low energies, there are two distinct peaks 
in the experimental data: the quasi-elastic peak at low 
photon energy and the pion-production peak at higher 
photon energy. At lower energies, we overestimate the 
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FIG. 10 Same as fig. [9] at an incoming electron energy of 
1080, 1200 and 1500 MeV. 



experimental data in the quasi-elastic peak and give a 
good description of the pion production region. Note that 
especially at low photon energies, nucleon-nucleon corre- 
lations and particle-hole excitations (see e.g. (j4pl : l4ll )) 
and the broadening of the A resonance are expected to 
play an important role. Up to now, we neglected these 
effects in our model which might explain this discrep- 
ancy. Going to higher energies we observe a better agree- 
ment in the quasi-elastic peak and a lack of strength 
above the one-pion production peak. Here, we are at 
the present stage simply missing more energetic channels 
as e.g. j* N — > Nttt:. The overall agreement improves 
when the nucleon potential is taken into account. 
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V. SUMMARY 

We have studied pionic double charge exchange on dif- 
ferent nuclear targets ( le O, 40 Ca and 208 Pb) in the A 
region (£km = 120, 150, 180 MeV) with a semi-classical 
couple channel transport model (GiBUU). Furthermore, 
we compared the results of our model with the exten- 
sive set of data taken at LAMPF (f2~H ). achieving a good 
agreement, not only for the total cross section, but also 
for angular distributions. The scaling of the total cross 
sections pointed out in {H) could be reproduced. How- 
ever, we found that two important effects that break this 
scaling: neutron skins and Coulomb forces compensate 
each other. We have shown in section IIII.BI that the 
DCX cross section is very sensitive to the size of the 
neutron skin. A precise measurement of DCX at for- 
ward angles combined with a theoretical analysis could 
bring quantitative results on the neutron skins of heavy 
nuclei. We conclude that the implementation of pion re- 
scattering and absorption in the GiBUU transport model 
successfully passes the demanding test of describing dou- 
ble charge exchange reactions. Thus the semi-classical 
approach is well suited to describe pion dynamics in nu- 
clei for pion kinetic energies greater E^m ~ 30 MeV. 

As a second application we have presented our treat- 
ment of electron scattering off nuclei. In this problem we 
have employed in-medium kinematics and parametrized 
the form factors for QE-scattering according to Budd et 
al. H2) and Krutov et al. (f33h and used the 7r-production 
form-factors of the MAID([36|; [U l38l ) analysis. We have 
evaluated so far only inclusive cross sections and com- 
pared to data by Anghinolfi et al. ((391 ). In comparison to 
the data we found good agreement in the pion produc- 
tion energy regime, whereas the quasi-elastic scattering 
is overestimated for low photon energies. In future we 
plan to extent the description to exclusive channels, such 
as e.g. single 7r production, and plan to include more 
energetic channels like e.g. -k-k production and a more 
realistic spectral function for the nucleon |42|; EF 
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